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This paper presents a physics-based mathematical model for a fuel cell 
hybrid vehicle system. The performance design of fuel cell hybrid electric 
vehicles (FCHEVs) is an area of interest in transportation applications. 
FCHEV is a combination of a proton exchange membrane fuel cell 
(PEMFC), with a battery and associated DC/DC and DC/AC converters. 
Suitable batteries are used in FCHEV acts as a backup system with efficient 
energy management. The battery is designed for fast power transfer during 
transient response and constant performance without hydrogen. The power 
management schemes discussed are consistent with the most common 
modern power management techniques used in fuel cell vehicle applications, 
including: the state engine control strategies, rule-based fuzzy logic 
strategies, classical proportional integral control strategies, frequency fuzzy 
logic decoupling/control minimization strategies. Therefore, we tested and 
evaluated the reliability of the model using the MATLAB/Simulink 
environment, with mixed results. 
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NOMENCLATURE 
T, : The total traction force Vactivatic: Voltage loss due to the reaction rate 
at electrode surface 
Tyf : The traction force fronts the wheels y? : Open cell voltage of the fuel cell 
T : Traction force of the rear wheels F : Faraday constant 
Tm : Motor torque To : The standard temperature (25 °C) 
Ra : Force has the air resistance on the vehicle AS? : Entropy change (considered 
constant = -164 J/mol. K) 
fr : The rolling resistance coefficient R : Universal gas constant = 8.31 
J/mol. K 
Re : Total strength of the resistors Pnz» Po, : The pressure of hydrogen and 
oxygen 
Rey : Total resistance force of the Wheels D : The air pressure on the electrode 
on the ground surface 
Pn : The engine required to drive the vehicle Po : Atmospheric pressure 
at a speed 
Py : Strength of resistance to the slope Pair-com: Power needed to compress air 
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P, : Total power Pair-cir : Power consumption of circulating 
air pump 
Net : Transmission system performance Mair  : Mass flow 
M, : Global mass of the vehicle m : Constants in the mass transfer 
voltage 
ig : Gear ratio of transmission Ne : Cell number 
ig : The gear ratio of the final drive Irc : Output current density [A/cm2] 
Tj : The effective radius of the tire I, : Internal current density 
Pa : Mass volume of air lo : Exchange current density 
Nt : The efficiency of driveline from the power Is : Output current fuel cell 
plant to the drive wheels 
Ar : The front surface of the vehicle Lax : Maximum current fuel cell 
Cp : Aerodynamic drag coefficient T : Internal resistance fuel cell 
Vo : Basic speed of the electric motor "pump : Efficiency of the air pump plus 
motor drive 
Vy : Wind speed and direction of movement Ny : Efficiency voltage 
V : Vehicle speed NF : Faradic efficiency 
N : Motor speed in revolutions per minute Nsys : Efficiency system 
AG : Gibbs free energy Tota; : Total efficiency of the fuel cell 
Gn : Free energies in species “n” of products “a? Ufy, : Use of hydrogen fuel cell 
Gm : Free energies in species “m” of product AH : The reaction enthalpy at the 
reactant absolute temperature T 
Z : Number of moles of electrons exchanged FCHE : Fuel cell hybrid electric vehicles 
for the formation of one mole of water V 
B : Charge transfer coefficient ESS : Electrical storage systems 
y : The ratio of specific heats of air (= 1.4) PEMFC : Proton exchange membrane fuel cell 
Vath : The theoretical nominal voltage AFC  .: Alkaline fuel cell 
Veoncentration : Voltage loss due to reduced gas PAFC : Phosphoric acid fuel cell 


concentration or delivery of large 
quantities of oxygen and hydrogen 

V : Output voltage of fuel cell SOFC : Solid oxide fuel cell 

Vonm : Ohmic voltage drop due to resistance of MCFC : Molten carbonate fuel cell 
proton flow in electrolyte 


1. INTRODUCTION 

In the ambition to deal with environmental pollution and the limited supply of fossil fuels, 
automotive industries have embraced the objectives outlined by international and regional agreements for 
environmental protection and pollution control and have shown increasing interest in fuel cell hybrid electric 
vehicles (FCHEVs). This orientation is evidenced by the recent agreement of the European Parliament 
adopted on Wednesday, June 8, 2022, which aims for a ban on combustion-powered cars and vans from 2035 
within the European Union. Rapid advances in electronics technologies and the great potentials offered by 
fuel cell have enabled a developed future for FCHEVs [1], [2]. 

The potential of FCHEVs remains enormous in terms of cost reduction and improved durability. 
The implementation of a second energy storage participates to decrease significantly carbon dioxide 
emissions in the vehicle [3], [4]. For that reason, the hybridization process plays a major function in the 
rational consumption of hydrogen and reducing the greenhouse gases [5]. 

Typically, FCHEV is composed of two different energy sources in perspective to use their 
advantages to enhance the performance of the FCHEV, the storage system has proven its reliability and 
technical feasibility. Compared to conventional combustion engines, they provide efficient electricity, lower 
noise and virtually zero emissions. That’s why it seems necessary to make a combination between fuel cells 
and other energy storage component such as lithium-ion batteries or super-capacitors. This blending can 
optimize the exchange of energy and extending the autonomy of the fuel cell systems by taking advantages 
from the battery charge. The traditional way to evaluate hybrid system performance is to create a simulated 
model system of the vehicle, which allows the test of many different driving cycles, the tire dynamics and 
breaking effects on vehicle performance haven’t generally been considered. 
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2. METHOD 

As represented in the Figure 1, the system we will use in this paper contains 4 global blocks or 4 
main functions. The driver block represents the standard driving cycles in order to simplify the approach. 
Concerning the control system, it aims to establish an energy management system that can be adapted to the 
driving conditions and which will then be transmitted to the vehicle system’s block for application. The 
speed and the energy profiles issued from the vehicle system are viewed at the result’s block. 


DRIVER CONTROL VEHICLE 
SYSTEM : SYSTEM SESULIH 
(DRIVE CYCLE) 


Figure 1. Representation of the fuel cell hybrid vehicle 


2.1. Control system 

The strategic management of Energy tries to explain how energy is distributed between batteries and 
fuel cells. They can optimize the fuel cell capacities and consumption depending of energy demands. When 
developing and evaluating energy management strategies, it is necessary to understand the ultimate limit to 
improving that fuel cell stack aging in relation with driving manners. Jung ef al. [6] are proposed an energy 
management strategy, based on maintaining efficient fuel cell stack and cell behavior by using a hybrid 
controller. Cell and demand power management between fuel cells helps to guarantee a transitory behavior of 
the battery pack with more efficiency and a slow aging [7]. They developed a power distribution that could 
account for the loss of recovered energy loads. 

The optimization process declined in many developed algorithms is aimed principally on reducing 
the power consumption of the propulsion system. In spite of that, in the future, regenerative braking energy 
of the electric motor must be included simultaneously with the FCS output power for recharging the 
battery [8]. Therefore, when considering regenerative braking energy, the power consumption of the car must 
take into consideration the battery consumption level. Beyond this particular problem, a control strategy 
developed specifically for FCHV is presented in [9]- [14]. 


2.2. Vehicle system 

The basic principles of vehicle operation represent vehicle behavior mathematically as an interactive 
mechanical system. A vehicle in a complex ecosystem comprising a number of interconnected components. 
Complex mechanical and mathematical knowledge is required to adequately describe its behavior. As shown 
in Figure 2, a representation of the vehicle's resistance including tire's ones to its motion. 


Figure 2. Vehicle body longitudinal dynamics 


The formula for the mechanical power of the engine requisite to impel the vehicle at speed V: 


VTt 


nt 


P, = (1) 


with T, : the total traction force T, = Tyf + Tyr. The dynamic equation of vehicle motion is expressed by 
Newton's second law: 


dV _ Tt-Rt 
dt 6M, 


Q) 
with 6 : the factor of the mass. 
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The resistances of the vehicle that oppose its motion include the rolling resistance of the tires. The 
expressions of F, of air resistance force on the vehicle, with F;,, the global resistance force of the wheels on 
the ground (Rt; = Rray + Rrar) are the following: 


_ PaAfCp (VV)? sgn(V+Vy) 


Ra = eo 3) 

Re = M, g f, cosa (4) 
when the vehicle moves up (down) a hill, a force of resistance to the slope caused by its weight: 

P, = M, g sina (5) 
the total force of the resistances to movement is: 

Ri = Ra + Ry + Py (6) 


consequently, we neglect V, and replace in the dynamic equation of motion of the vehicle [15]. 


dv Tt-5paAgCpV? +My g sina+My g fy cosa 7 
dt 5My V) 


So, if a is small (c.-à-d. cos a, sina = tana ): 
dV | 1 2 
T, = oM, + 5 PaAf Cp V + M, g (f. +k) (8) 


with k = tana. This gives the total power that requires the vehicle to accelerate from zero to that speed V in 
second: 


T (V? + v2) + = pa Ar Cp V? + M, g (f, + k)) (9) 


1/2 
P, = + (2M, gVf-+ 


2.3. Transmission 

As modeled in Figure 3, the transmission and chassis simplified system is presented. Electric cars 
use a single gear, which makes acceleration and braking easier. Wheel traction (T+) and vehicle speed (V) can 
be expressed using the following equations: 


Tmigiont 
T, = ERU 


"s (10) 
.nNrgq 
~ 30igio (11) 


Cà Control. Brakes . i (l&——À»1 
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| | 
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Figure 3. Model of the "Transmission and Chassis" system in Simulink SimPowerSystems 
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2.4. Motor 

The output voltage of the suggested transformer is sent through an inverter to an electric motor for 
vehicle propulsion. Electric vehicles play an important role in FCHEV. The right electric vehicle can 
significantly reduce the size of fuel cells and its costs consequently. In the beginning, most automakers resort 
DC motors for electric vehicles. The brushes and rotating devices make DC motors more expensive to 
maintain and inefficient operating [16]. Currently, the technology using motors BLDC with a permanent 
magnet is the most used in FCHEV applications, like the model presented in Figure 4, due to their ease of 


control, high reliability and robustness. 
m m.i 


m. w @4—4——-v 


Torque Sensor 


Motor 


Figure 4. Model of the electric motor in Simulink SimPowerSystems 


2.5. Fuel cell 

Fuel cells are supposed to be low-emission generators for stationary and mobile applications. There 
are different types of fuel cells using different electrolytes, the operating temperature, and most importantly, 
the field of application. They can be divided into four categories, details in Table 1: solid oxide fuel cells 
(SOFC), (MCFC) molten carbonate fuel cells, phosphoric acid fuel cells (PAFC) and (PEMFC) polymer 
electrolyte membrane fuel cells. All these systems, PEMFC dominates the automotive industry with its 
minimized operating temperature and fast start-up time [17]. 

PEMFC use a perfluoro sulfonic acid membrane (usually Nafion) as an electrolyte; their 
performance is highly dependent on their hydration state [18], [19]. They are the famous fuel cells integrated 
in automotive applications. This is greatly because they work at the low temperature so by a quick start the 
full load can be achieved easily. 


AG = È product Gn = È reactants Gm (12) 


Where G, and Gn are the free energies of product species n and reactant species m. 

As shown in the Figure 5, ion exchange membrane fuel cells use a reaction between hydrogen (or a 
hydrogen-rich fuel) and oxygen (or air) to create an electricity energy based on two redox reduction 
reactions: 


anode H, > 2H* +2e7 


1 
cathode 2H* + 502 +2e”- > H,0 
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Table 1. Table representing several fuel cells 


Type of FC Electrolyte Operating Chemical reaction 
temperature 
Polymer electrolyte Solid organic poly 60-100 Anode: H, > 2H* + 2e7 
membrane fuel cells perfluoro sulfonic acid ; 1 " _ 
(PEMFC) Cathode: 202 +2H* *2e€ > H0 
S 
Cell : H, +50, > H0 
Phosphoric acid fuel Liquid phosphoric acid 175-200 Anode: H,  2H* + 2e7 
lls (PAF aked i tri 1 
SEIS (PARC) SORS ER RU de Cathode: 0, + 2H* + 2e~ > H,0 
a e P Ó—— 
Cell : H, +50, > H,0 
Molten carbonate fuel Liquid solution of Li, Na, 600-1000 Anode: H, + CO, > H,0 + CO; + 2e7 
lls (MCF K carbonat ked i 1 
cells (MCFC) carbonates soaked in a Cathode: = 0, + CO, + 2e~ > CO, 
matrix 2 
Cell : H, t 20; + CO; > H0 + CO; 
2 
Solid oxide fuel cells Solid zirconium oxide 600-1000 Anode: H, + 0; > H,0 -2e^ 
(SOFC) with small amount of . 1 E _ 
yttria added Cathode: 29: +2e 0; 
—— ammm 
Cell : H, +302 > H,0 
Charge 


Entrée de 


Entrée de 
lons de 
combustible RT Eydrciséns p nt 
H+ 
02 

Sortie de Sortie de 

combustible carburante restant 
restant 


et gaz produits 


électrode Anode Electrolyte électrode Cathode 


Figure 5. The reaction in a fuel cell with an exchange membrane 


The overall chemical reaction is described by the equation [20]: 
1 
H, + 2 0,  H;0 


This reaction is spontaneous (it does not require energy to do so). The energy delivered in reaction to form 
water, called the Gibbs free energy (AG), conditioning by temperature, the operating pressure of the cell, and 
the state of the water at the outlet (which may be in liquid or gas form) [21]. Therefore, the theoretical 
nominal voltage was established from the Nernst equation: 
— Ae = 0 APER pd 

Vatn = — 2 =V + (T +T (—) i (13) 

The output voltage of the PEMFC Vs is defined as an outcome combination of the cell current, the 
pressures of the reactants, the ambient temperature of the fuel cell and the humidity of the membrane, 
represented by these equations: 


V, = Vatn — Vactivation ^ Vonm — Veoncentration (14) 
with: 
RT Ifctln 
Vactivation x Ne 2BF m (15) 
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In fact, fuel cells require auxiliary materials to function. Auxiliary equipment is mainly composed of 
air circulation pump, coolant circulation, pump, fan, oil pump and the electrical control component. In 
auxiliary units, airflow pumps are the biggest energy consumers’ power consumption of the circulating air 
pump including its drive motor can reach 10% of the total fuel cell output. 

Auxiliary equipment is considered less energy consumers comparing to the circulating air 
pumps [22]. In fuel cells, the gas pressure p at electrode surface is often above atmospheric pressure pọ to 
reduce the voltage drop. From thermodynamics, the power required to upgrade the air pressure from 
compress po to another pressure p value with Mair mass flow can be presented as (18): 


yi 
Pair-comp = hare R. T © di. 1) (18) 


0 


with T: the compressor inlet temperature in Kelvin. 
In the calculation of the power consumption of the circulating air pump, the total power consumed is 
presented as (taken into consideration the energy losses in the air pump and the motor drive): 


Pair-com 
Pair—cir ERE (19) 
Npump 
— Fuel cell efficiency 
The voltage efficiency, is the ratio between the real voltage of the fuel cell and Vatn battery (given 
by the Nernst equation): 


Vs 


qo = y —— 


Vnth 


(20) 


The faradic efficiency np, is the ratio between the output current of the CAP and the maximum current: 


I. 
ns = —— 


Imax 


(21) 


the system efficiency fsys, is the ratio of the net power delivered by the fuel cell generator (the system) to the 
power generated by the fuel cell: 


Neys = l4 (22) 


Pgross 
Total system efficiency Ntotale, given by the (23): 


ZNsys F Uf uS Vs 


Ttotal — —— >, (23) 


AH 


with large currents cause inefficiencies due to large stack voltage drops, while very small currents cause 
inefficiencies due to the increased percentage of power dissipation in auxiliary power dissipation [23]. 

A model of the fuel cell system was created in MATLAB/Simulink using the SimPowerSystems as 
presented in the Figure 6. The PEMFC stack produces a low unregulated DC output voltage. Traditional 
boost converters are used as the power electronic interface, and boost converters may not be compatible with 
high power applications due to their low current and thermal management systems [24]. To deal with these 
difficulties, various high voltage gain DC/DC converters have been adopted [25], [26]. 


2.6. Battery 

Electric vehicles have the advantage of being quiet, extremely low or even zero-emission. 
Nowadays, battery technology is the principal serious issue in electric vehicle development. Vehicles with 
internal combustion engines. To enhance a fast and a sure development of electric vehicles, one effective way 
is to get effective solutions for battery issues, and also to solve the battery problem. 
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Today, fuel cell hybrid electric vehicles (FCHEVs) are taking more importance in the automotive 
industry. Generally, in a fuel cell hybrid vehicle, a battery is made for use in parallel to the fuel cell. The 
battery is responsible for starting, cold cranking, electrical braking surges and energy storage ministry of 
electricity [27]. 

The batteries, studies and chosen in this paper are lithium-ion batteries because of their higher 
energy density and better efficiency than other battery types such as lead-acid, nickel-cadmium, or nickel- 
metal hydride [28]-[30]. This presents certainly a great potential for automotive or aerospace applications. 
Lithium-ion batteries are electrochemical energy stores and the source of electrical energy. Lithium-ion 
battery cells can be presented using the equal circuit according to [31]. This model of battery is characterized 
by sensible accuracy, simple parameter determination and short service life. 


Boost Converter Control 


Current Sensor 


Fuel Cell 


Vehicla . Vv @— 


Control . iCmd 
Control. TFuelCellCmd  —— »| Temi. 


Figure 6. The fuel cell system created in Simulink using the SimPowerSystems 


3. RESULTS AND DISCUSSION 

To highlight the performance of the suggested model of the hybrid vehicle, Figures 7 and 8 show the 
results of the simulation in the NEDC and the FTP driver cycle: These figures show the hydrogen 
consumption every 100 km, the power is communal between the battery and the fuel cell in the step of 
storing negative power recovery in the battery. Depending on the load of the driving cycle described in the 
Figures 7 and 8. The energy losses reach a global optimum based on the reference power profile. 

The simulation results show that after the peak hydrogen consumption reaches 24 L/100 km in the 
starting phase of the vehicle, the hydrogen consumption drops significantly and gradually stabilizes at the 
value of 2.5 L/100 km, as show in Figure 8 which is very meaningful. Driving cycles aims to normalize the 
simulation and the illustration. Additionally, it appears clear that the impact of tire wear on powertrain 
performance and energy efficiency cannot be ignored during severe acceleration and braking. Also, vehicle 
control objectives were identified in maximizing the enormous benefits associated with the functional 
integration either of fuel cell and battery technologies. 
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Figure 7. The results of FTP driver cycle 
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Figure 8. The results of NEDC driver cycle 
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4. CONCLUSION 

Detailed mathematical and physical models of battery-equipped FCHVs have been presented in this 
paper. This model is suitable for studying the behavior and the dynamic performance of the propulsion 
system during heavy acceleration and regenerative braking. It contributes to a better understanding of the 
achievable efficiency of the vehicle in the "real world". Powertrain transient and steady-state losses are 
quantified in predefined test cycles containing higher accelerations, travel and braking values than 
conventional legal driving cycles. Simulation results show that the impact of losses on powertrain 
performance and fuel efficiency cannot be ignored during rapid acceleration and braking. 
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